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Abstract 

The theoretical computing of special values assumed by the hypergeometric functions has a high inter- 
est not only on its own, but also in sight of the remarkable implications to both pure Mathematics and 
Mathematical Physics. Accordingly, in this paper we continue the path of research started in and [20] 
providing some contribution to such functions computability inside and outside their disk of convergence. 
This is accomplished through two different approaches. The first is to provide some new results in the spirit 
of theorem 3.1 of [20], obtaining formulae for multivariable hypergeometric functions by generalizing a well 
known Kummer's identity [IS] to the hypergeometric functions of two or more variable like those of Appell 
and Lauricella F^^ 

In the second part, using some reduction schemes of hyperelliptic integrals due to Goursat [5] we evaluate 
Appell and Lauricella F^' functions in some particular occurrences and in their analytic continuation. 

Keyword: Reduction of Ifyperelliptic Integrals; Complete Elliptic Integral of first kind; Hypergeometric 
Function; Appell Function; Lauricella Function 

1 Introduction and aim of the paper 

The theoretical evaluation of the hypergeometric functions has became of great consideration not only in order 
to provide a set of exact values for testing the numerical algorithms of their computation. But the best chiefly 
comes from being such aim accomplished by means of mathematical technics of remarkable power and interest. 
In fact it was by studying the hypergeometric series that two of the most important chapters of the functions 
theory were originated at the end of the XIX century, namely that on the linear differential equations (I. L. 
Fuchs), and the automorphic functions (H. Poincare, F. Klein). Furthermore F. Klein has shown that from the 
hypergeometric equation just solved by the hypergeometric function, really comes out the majority of the linear 
differential equations of Mathematical Physics, if not all. 

Just the uncountable applications of such functions to the Mathematical Physics have attired the numerical 
analysts in order to improve the quality, reliability and celerity of their effective computation. Coming to more 
recent developments, we can cite e.g. three papers by Joyce and Zucker, [T31 [23 [T3], where they employ for 
computational purpose the theory of singular elliptic moduli. 

In order to provide some contribution to their computability, almost as a reference for some values of them 
inside and outside their disk of convergence, in this paper we continue the piece of research started in (19j and 
[20] following two aims. The first is to provide some new results in the spirit of theorem 3.1 of [20] where we 
proved that, if 2& > a > 
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What above follows from the double evaluation method (namely in terms of Gamma function and using the 
integral representation theorem, see (|2.2p below) applied to the improper integral 
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Similarly, here we obtain formulae involving hypergeometric functions of two or more variables (Appell Fi and 
Lauricella F^-* ) using two sequences of integrals 

AJa,b)= / -rdt, BJa,b)= / -rdt. 

In this way we succeed in generalizing the Rummer's identity (see [TSl, [3] pp. 9-10, or Cor. 3.1.2 p. 126 in [T]): 

\ r(l + a-6)r(l + f) 




r(l + a)r(l + f -6) 



(1.1) 



concerning the Gauss function, to hypergeometric functions of two or more variable like those of Appell and 
Lauricella F^-* . 

In the second part, using some reduction of a class of hyperelliptic integrals due to Goursat [5] we evaluate 
Appell and Lauricella F^'' in some particular occurrences and in their analytic continuation. This continues 
our previous research of [12] and [5D] where we obtained several identities connecting tt, elliptic integrals and 
hypergeometric function using some reduction schemes of hyperelliptic integrals due to Jacobi |12] and to 
Hermite [11]. 



2 Notations 

In this paper whenever 2F1 denotes a well-known hypergeometric series: 

(«)™ (bin ^" 
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(2.1) 



it also denotes its analytic continuation on C \ [1, 00) via its integral representation theorem Re(a) > 0, Re(c - 
a) > 0, |x| < 1: 
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In (|2.1I) {a)m is a Pochhammer symbol: 



r(a + m) 
T{a) 

We also consider the multivariable extensions of 2F1: the Appell Fi two- variables hypergeometric series, see [5], 
defined for |a;i| < 1, \x2\ < 1: 
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the analytic continuation for the Appell's function on C\ [1, 00) x C\ [1, 00) comes from its integral representation 
theorem: if Re(a) > 0, Re(c — a) > 0: 
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a; 61, 62 
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Its n- variable extension is the Lauricella hypergeometric series, introduced in |16) : 



■ du. 



a; 61, ... , 6„ 



00 00 



Xi , . . . , X7J 



(a),„i+...+™„(&i)„i •••(6„) 



mi— nin—Q 



(c)„,.i+...+„j„mi! • • -TTln! ^ 



whose integral representation theorem for Re(a) > 0, Re(c — a) > allows its analytic continuation on the 
n-fold cartesian product of C \ [1, cx))" 
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We will also use the order reduction formula for Lauricella functions 



a;bi, . 
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whose proof was given in |20| Lemma 1.1 therein. 
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3 Special values from Eulerian integrals 
3.1 Evaluation in the boundary of the unit disk 

We will present here a generalization of the Kummer identity to both Appell and Lauricella functions. Our start- 
ing point is an elementary proof of we discovered, at least we presume, throughout a two-fold evaluation 
of the definite integral: 

Ma.b)^f^^^At (3.1) 

where, for the sake of convergence of (jS.ip . we assume Re(a) > and Re(6) < 1. First we change variable in 
p.ip putting = u; it follows that 

ir(f)r(i-6) 



On the other side, if we write p. ID as 
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From the integral representation theorem, equation (|2.2p . we get: 

, , r(a)r(l -b) ( a;b 

T(l-a-b) \ l + a-b 

By comparing both expressions for ^2(0, b) and recalling that 

r(§) _r(i + §) 
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(3.2) 



(3.3) 



2r(a) r(i + a) 

we infer (|l.ip which was originally given by Kummer in |15) . 

To obtain for the Appell Fi two variables hypergeometric function something analogous to formula (11.11) . 
we consider the integral, where we changed the exponent 2 in p.l[) to 3: 
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Again putting = u we evaluate p. 41) through the Eulerian integral of the first kind as 



^3(0,6) 



i r(f)r(i-b) 
3r(i + t-&) ■ 



On the other side we can factorize 1 ~ obtaining 



Aaia.b) 



dt^W^F, 



lo {l-ei''H)\l-e-i''H)'' r(l + a - 6) 
Comparing the expression for A3 (a, b) we get 
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Notice that we use the trivial identity, n G 



nr(a) r(l + a) 



taking n = 3. 

Continuing, from the integral 



1 ^a-l 
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factorizing \ — t'^ we get the Kummer-like formula for the Lauricella F 
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We see that the general relation for n e N, rt > 2 is 

A„(a, 6) = 



- 1, i, -2 = 
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Thus we can state our result in general. In order to simplify the notation whenever in a Lauricella function the 
n parameters as b are all equal, 6i = • • • = 6„ = 6, we put 
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a; 6i, . . . , 6„ 



Xi , . . . , 
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the repetitions number of b is denoted by the apex of the Lauricella function. 

Theorem 3.1. Let n he a positive integer > 2. For k ~ 1, . . . ,n — 1 let uj^^n ~ so that wj!^ — 1 and 

^k.n 7^ 1 for = 1. Let < Re{b) < 1 anrf Re{a) > t/ien 



.(n-l) 



a, b 
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Proof. The statement follows from the integral representation theorem (|2.3p from which we can establish relation 
p.lOp below that expresses integrals A„(a, b) in terms of ri — 1 variables Lauricella function F 
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Thesis follows comparing p.lOl) with 



l + a-b 



(3.10) 



□ 



Remark 3.2. In the case n = 2 when a = 1 and & = 1/2 we have the elementary value ^2(1, 1/2) — n/2 which 
drives to 

1,1/2 



2F1 

according to equation (14) section 2.8 p. 102 of [7] 
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3.2 Values outside the unit disk 

Consider the sequence of integrals 

B„{a,b) 



-dt 



(3.11) 



lo (1 + i")' 

where we assume a>0,6>0, na>6to ensure convergence. The integral at the right hand side of p.lip can 
be computed by means of the Mellin transform in fS^ p. 310 equation 21, see also [TU] entry 3.194-3 p. 313, so 
that 

On the other side, we can also evaluate i?„(a, b) using the integral representation theorem (|2.3|) . In fact change 
of variable t = {1 — u)/u drives to 



Bnia,b) 



1 y^7ib — a—l 



(1 - uY 



(u" + (1 - m)") 



-du. 



(3.12) 



Observe that for n even the polynomial u" + (1 — u)" has degree ri, while for n odd has degree n — 1. In any 
case the roots are distinct and the factorization holds: 



<i-«r-n(i-i") 



where ak is root of u" + (1 — u)" = 0. These roots are evaluated, through De Moivre formula, as follows. If n 
is even, say n = 2m, we have 
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If n is odd, say n = 2m — 1 



(2fc-l)7r 
2m-l 



2fi + cos(^^) 
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, fc = 1, . . . , 2m — 1, fc 7^ m. 



We are in position to state our second result. 
Theorem 3.3. Assume that a > 0, b > 0, nb > a. If n = 2m is even then 
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where, from p.lSp for fc = 1, . . . , 2m.' 
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If n ~ 2m — 1 is odd then 

p(2m-2) / (2"'' - 1)^ - a; 

° \ (2m -1)6 
where from (I3.17P for k — 1, . . . , 2m — 1, k ^ m 
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Remark 3.4. In the even case n = 2m it is possible, using Lemma 1.1 of [20 , reduce the order of Lauriceha 
function appearing in p.lSp . obtaining: 
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with arguments given for k = 1 , . . . , 2m — 1 by 
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Thus for TO = 1 we obtain, as a particular case, theorem 3.1 of |20| . 

Proof. Using integral representation theorem, formula (12. 3p in equation p.l2p we get, if n = 2to is even 
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where we denote with Xi, . . . ,Xn the reciprocal of the n roots of equation it" + (l — m)" = given by (j3.16p . 
While if n = 2to — 1 is odd we have 



T{nb) 
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but here root x^+i = Xm is skipped. Theses p.lSp and p.l7p follow comparing respectively p.l8p and p.l9p 
with (|XTT|) . ' □ 



4 Some special cases 



When in some special situation the Eulerian integrals yl„(a, 6) and i?„(a, 6) can be evaluated in terms of 
complete elliptic integral of first and second kind we are able to establish, from our formulae p.9p and p.l7p , 
new relationships between hypergeometric functions of one or more variable at some special values and certain 
complete elliptic integrals. This relationships arise for singular moduli of the elliptic integrals, according to the 
relevant theory studied by several authors, see [331 IMl HI IS]- We will use some integrals taken from [TO] and 
some less known integral studied in the past by some eminent mathematicians of XIX century, Legendre |17] 
and ^ISj, Richelot ^1^, Serret [53], integrals, as far as we know, almost forgotten by the recent literature. 

Starting with the most elementary elliptic integrals, notice that in any occurrence we also express the integral 
under investigation in terms of our families A„(a, b) or -B„(a, b) 
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Using (Pl^ and ([51^ we see that 
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Using Lemma 1.1 of [50] in (14. 3p we also establish values of Y)j when one of the variables is in the positive real 
axis in analogy with our results for the Gauss 2F1 exposed in [20j theorem 3.2 equations (27) to (30) therein 
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Now we turn our attention to some integrals from classical repertories TO' and (61. We begin with integral 
3.183.2 p. 313 of ilOj 



so that comparing with (j3.3p we get the formula 
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From entry 3.184.1 p. 314 of [TO] 
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we obtain 
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Using formula 3.185.2 p. 314 of [TO] 
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We end the examples taken from [TO] with entry 3.185.4 p. 314 
which drives to 
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Consider entries 576.00 p. 256 and 578.00 p. 258 of ^6] 
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Comparing with (IXTO)) and with (|XT5)) we obtain 
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and 
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Observe that employing p. 191) in (I4.12p we evaluate F^-* when one of its arguments is 2. Putting 
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Now we consider some integrals studied by Legendre, Richelot and Serret. Legendre, [TB] P- 383, proved 
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thus comparing (|4.14p with p.lOp we get 



Richelot in 



D 



evaluated 



1 + i . 1-i 1 + i . I - i 



V2 ' ' V2 



^8(3, i)= /' 
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thus comparing (|4.16p with p. 101) we obtain 



3; 



. 1 



1 + i 1-i 1 + i 1 - i 
h ^,-1, 



V2' ' V2 



V2 ' ' V2 



16 



2V2 



K ( V2- 1 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



To establish further identities, we recall two huge formulae due to Legendre, that, as far as we know, are quite 
unrecognized by the mathematical community of our age. In his famous Traite, [18] p. 377, equations (z) 
therein, he established two remarkable formulae, the first, for 2a < n 
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the second, for n/2 < a < n, is 
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: dx = 



2n 



n(2a — n) sin(— a) 
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We do not want to report here the path followed by Legendre to establish (|4.18p and (|4.19p . a very interesting 
topic, with the original source now easily accessible. We highlight that (|4.18p plays a role in order to evaluate 
the integral (I4.14p so that it can be used also for evaluating the integral (j4.16l) . 

Now by formulae (j4.18l) and (I4.19P let us pass to some evaluation of elliptic integrals related to (j4.14l) and 
(OBI) : 
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Thus taking advantage of p.lSp in (|4.20p and (I4.2ip we see that 
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with xi, . . . , Is reciprocal of the roots of + (1 — u)^ = as given by p.l8p . namely: 
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We can also reduce the order of the Lauricella functions in (j4.24p and (|4.25p using Lemma 1.1 of [20] getting 

(7) 

special values of F)j with some of the arguments lying in the positive real axis 
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where 
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1 + I . 1 - i 



2:4 = 2, Z5 = 1 



1_M 
V2 ' 



,(7) / ^' 2 



15 



11 
2 



15 
2 



72'"' 72 ' ^' 
1 + i 1 - i 



1 + i 1 - i 



SIStt 



1 - i 

71" 



%/2 ' ' V2 y 1024%/2K (V2 - 1) 
1001 (2 + ^/2) TT 



1 + i . 1 - i 



V2 ' ' V2 ' ' V2 ' ' V2 / 16384K(V2-1) 



(4.28) 
(4.29) 



We add to our list of integrals one more, considered by Serret, [23j p. 65 and Legendre, ^/Tj chapter XXX 
pp. 201-202, but (again) ignored by classical handbooks: 
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Eventually from p.l9p we get an evaluation in the analytic continuation of fI^'' . Putting 
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In Chapter XXX section 164 pp. 205-206 of [TJ Legendre evaluated an integral similar to (j4.30l) 
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So using p.lOp we get the identity 
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Furthermore in this special situation it is possible to use equation (|2.4p . i.e. Lemma 1.1 of [5D], also for an 



integral of the A family, in (j4.34p obtaining 



1 %/3 . 



5 Special values from reductions of hyperelliptic integrals 

In a 1885 article [5], Goursat, using the change of variables: 

t3 + at + b 



(4.35) 



obtained the reduction formula: 



dx 



y/x (4(3a: - a)3 - 27{b + px)^) 

-ap-ib-2,t^{p-2t) 

^{at + h + t^) [ap + 36 + 3i2(p - 2t))^ (4(ap + 36) + 3pt2 + 3^3) 



dt 



(5.1) 



Once a suitable interval of integration is chosen, it is possible to remove the sign uncertainty and simplify (|5.ip 
obtaining 

(5.2) 

y^x (4(3x-a)3-27(6 + pa;)2) 7 ^(at + 6 + t3)(4(ap + 36) + 3pt^ + 3*3) 

If we make the simple choice to take a = p — and to integrate in [(6^/4)^/3^ the left hand side of ([57 
we get 

da; dt 

' (5.3) 



dt 



v/27a;(4x3 - 62) 7(|)V3 ^3(t3 + 6)(t3 + 46) 

Normalizing the interval of integration we see that the most natural choice 6 = 2 is not restrictive, so we can 
use, as starting point to find our next identity, the integrals equality: 

r ^ -6 r dt 

Ji ^x{x3 - 1) " J I ^{t^ + 2){t^ + 8) 

The hyperelliptic integral at the right hand side of (|5.4p is a particular case of the family of integrals described 
by the following Lemma. 
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Lemma 5.1. Let 2n — m > 0, a, b > then 

dt 



"^(t" + a)(i" + b) 2n - m 



Fi 



2n-m . J_ J_ 

mn ' m ' 7n 

2n—m-\-mn 
mn 



(5.5) 



Proof. Condition 2n — m > ensures the convergence of the integral at the right hand side of (j5.5|) . Thesis 
foUows immediately from the integral representation theorem for the Appell Fi after the change of variables 
f=l/ui/". □ 



From (j5.4p and Lemma |5. II we infer the following new evaluation in the analytic continuation of Fi. 
Theorem 5.2. 



2.11 
3' 2' 2 

5 
3 



-2,-8 = 



1 



Proof. Using Lemma 15.11 in equation (|5.4p we see that 

dx 



But we also have 



V^i^' - 1) 



dx 



3Fi 



2.11 

3' 2' 2 

5 
3 



r3(i/3). 



dx 



r(i/3) 



/I v/^(x3 - 1) Jo 3 r(5/6) 

Thesis (I5.6p follows by equating (|5.7p and (15. 8p and using Euler reflection formula. 



(5.6) 

(5.7) 

(5.8) 
□ 



Remark 5.3. If, instead of using the Euler function to evaluate the integral in (j5.8p . we employ the elliptic 
integral of first kind, using jGj entry 244.00 p. 92 we get: 



Fi 



2.11 

3' 2' 2 



2,-8 



35/ 



' f^Lccos(2-V3),^^^' 



In order to express the Appell Fi not in its analytic continuation, but in the domain of convergence of the 
double power series which defines it, we present a slight modification of the Goursat reduction formula. 



Jo Vl - Jo 



6 



dt. 



(5.9) 



/o Jo ^it^ + 2){t3+S) 

Identity (j5.9p comes from the change of variable 

The hyperelliptic integral at the right hand side of (j5.9l) is computable using the Appell Fi hypergeometric 
function of two variables. In fact we have, in analogy with Lemma |5. II 

Lemma 5.4. Let a, b > then 

1-1 



dt 



V(*" + a)(*" + ^) 



ab 



Fi 



From Lemma 15.41 we get 



yit^TW+S) 4 ' 



1- J_ J_ 

i + i 

n 



1.11 
3 ' 2 ' 2 

4 
3 



1 1 

a' b 



1 1 

2'~8 



(5.10) 



(5.11) 



We can state a further identity involving tt which concerns a new evaluation of the Appell Fi. 
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Theorem 5.5. 



Fi 



1.11 

3 ' 2 ' 2 

4 
3 



1 1 

2'^8 



27^' 



(5.12) 



Proof. To get the thesis in this situation we evaluate the elhptic integral at the left hand side of (|5.9[) using the 
Euler Beta function: 

r^.^mm (.13) 

Jq Vl-a;3 3 r(5/6) 

Thesis (|5.12p follows from (|5.9p evaluating the integrals therein using (|5.1ip (right hand side) and (|5.13p (left 
hand side) and again invoking Euler reflection formula. □ 

Remark 5.6. If, instead of using the Euler Beta, we evaluate the elliptic integral at the left hand side of ([5 
using entry 244.00 p. 92 of jB] we get, instead of (j5.13p : 



Jo 7T=^ ^T3^ - j 



(5.14) 



leading to the elliptic analogue of (|5.12p : 



1.11 

3' 2' 2 

4 
3 



1 1 

2'^8 



= —= F arccos(2 - V3), 
3^ I 



\/6 + V2' 



(5.15) 



The family of integral identities pointed by the Goursat transformation is huge. Nevertheless the evaluation 
is intricate when some of the cubics in (|5.2p are complete. Taking for instance a = 0, 6 = p = 1 in (j5.2p we get 
the identity 



da; 



y^x{x - 1) (4x2 + 3a; + 1) Ji ^{t^ + 1) [t^ +f^+A) 



■.dt 



which can be rewritten as 



da; 



yj[l-x) (a;2+3x + 4) 7o y/W+T){W+t + T) 



3t 



■.dt. 



(5.16) 



(5.17) 



We need to employ the Lauricella F^'' in order to evaluate the hyperelliptic integral at right hand side of (j5.17p . 
To do it, first observe that: 



+ 1 = (1 + t) 1 - 



1 iVS^ 



1 iVl 



t I- 



1 iV?^ 



t = (1 + t)(l -cjt)(l -ujt) 



At' + t + l = {l + 2t)\l-\--^\t\\l-\- + ^\ t\ ^{l + 2t){l~et){l~-et) 



hence 



^(t3 + l)(4t3 + t + l) 



: di = - F 



(6) 



2- i 

^' 2 

3 



1, oj, —2, e, e 



(EHfc) 



We follow the same strategy to evaluate the integral at the left hand side of (|5.17p starting from the factorization 



3 



3 iV7 



+ix + A = A^-^-- ^ j xj \ l - y-- + ^ j xj = 4(1 - ax){l - ax 
leading to the integration formula: 

I' 

Jo ^(1 - a-) (a;2 + 3a; + 4) 
So we proved the following hypergeometric reduction formula 



dx 



Fi 



2 ' 2 

3 
2 



HSU) 
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Theorem 5.7. Defining: 
Then we have 







1 








2 








a, 


V 1 







(5.18) 



Remark 5.8. Of course, in evaluating the elhptic integral in (I5.16tip we can follow the elliptic way, |Bj entry 
243.00 p. 91: 



/' "^"^ = 7^ ^ f arccos (1(9- 4^2) ^ , \ J 8 + 5V2 

Jo J{l-x) (x2 + 3x + 4) 23/4 1^ V J J' 4^ 



lo v/(l - + 3x + 4) 23/4 

So, comparing (|5.17p and (|5.19p we get 



(5.19) 



2 



1, oj, —2, £, e 



3 



F Urccos ( y (9 - 4V2) j , i y^8 + 5\/2 j . (5.20) 




Remark 5.9. Observe that the Lauricella F^'' appearing in (j5.18p and (j5.20l) fulfills hypotheses of Lemma 1.1 
of [20^; thus in (|5.18p and (|5.20p the left hand side can be replaced with a lower order Lauricella. 

1 Z-iVl i + iy/l Z-iVl 3J-V7\ 
2' 4 ' 4 ' 4 ' 4 ) 

6 Conclusions 

The paper holds the proofs of a set of formulae providing the values in some special points -either inside or outside 
the unity disk- of some high order hypergeometric functions. Among them we can conclusively cite e. g. the 
following evaluations. First the p. 151) providing the Lauricella of order 2m; (j4.27p and the (14.28^ for Lauricella 
of index 7 through complete elliptic integral of first kind; the (j5.6p providing the Appell function through values 
of Gamma; equation (|5.15p provides an Appell function again through incomplete elliptic integrals of first kind; 
and so on. They have been theoretically obtained either by the double evaluation method or working on some 
established formulae for the reduction of hyperelliptic integrals. 
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